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Tit*  article  contain*  on  tviudyiln  of  the  proossaea  of 
mixing  In  a  turbulmt  at.iiW'idmre,  b.uivt  on  oystroAtiQ  ap¬ 
plication  of  the  nwthoda  of  Cm  theory  of  similitude,  Em¬ 
pirical  data  on  the  distribution  of  r»i rat  velouity  under 
various  oondtttona  of  temperature  atrauf toatton  ara  contr¬ 
ail  aod  anl  «  not  hod  la  nropoecd  for  computing  tha  auatauaoh  > 
oharaoumattoa  on  tha  haata  of  measuring  wt ni  voleolty  and  / 
temperature  | radiant,  y 

xmaDwnou 


Tha  questions  of  tho  phyaioa  of  tho  ground  layer  have  occupied  a  oonatdar- 
ablv  plaoa  in  aateorole steal  investigationa  during  the  poet  10-15  year a,  Tha 
laws  of  tha  proeeseea  lit  the  ground  layer  are  of  in  to  re  at  not  only  to  agro* 
meteorology,  vthioh  studio*  tha  offeot  of  a  "oetcorologioal  nodiua"  on  tho  growth 
of  vegetation,  but  they  alao  have  a  general  geophysical  significance,  since  tha 
dynamic  tnteraetlon  of  the  atroaphoro  and  tho  eubstrats,  the  "feeding*  of  tha 
atmosphere  by  aolatur#  and  heat,  la  realised  through  the  ground  layer. 

A  large  Mount  of  research  in  tho  field  of  ground-layor  phyalee  baa  boon 
done  at  tho  t'aln  Geophysical  Observatory;  the  works  of  S.  A.  Sapoahnikov*  [1], 

D.  L.  Lalkhtman  anl  A.  K,  Chuinovtkii,  (2],  U.  X,  Budyko  (3]  and  U.  T.  Timofeev 
{b]  are  well  taawn  to  Soviet  motcorologtets. 

Thte  research  hai  provided  valuable  observational  data  on  tho  dlatrlbutlan 
of  wind,  temperature  and  humidity  in  tho  ground  layer,  and  a  number  of  apooiilo 
propositions  have  bom  drown  up  on  tho  methodology  for  computing  turbulent  aua- 
tausch  characteristics  (Dudyko,  Laikhtman). 

In  thia  regard  Ciero  are  a till  a  number  of  debatable  quoations  in  tho 
thoory  of  grourv'-layor  mixing.  The  sin^Oort  system  of  the  "logarithaio  boundary 
layer,"  borrowed  from  technical  acrocfcmaaaics,  dasorlbea  quit#  well  tho  phenomena 
in  a  noutraUy  stratified  atmosphere,  and  la  supported  by  much  onplrical  data, 
however,  this  syotem  is  insufficient  for  describing  processes  In  a  real  ataos- 
phoro  where  the  tesperaturo  inhomogonoity  la  an  essential  factor  influencing  tho 
ilovelopmont  of  turbulenco.  Thia  latter  fact  (the  temperature  inhomogonoity)  dw- 
termlnos  tho  specific  nature  of  the  problem  of  atmospheric  turbulence  as  applied 
to  ground-layer  phyaioa. 

Tho  works  of  laikhtman  (5]  and  Budyko  [3),  as  well  as  those  of  a  number  of 
foreign  researchers  (Sverdrup,  R aaaby  and  liontgcoory;  see,  o.g.,  (6j)  have  boon 
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dovutrd  to  computing  th*  influence  of  tessera  turn  etratlflcaticn  on  turbulent 
exchango.  Th*  individual  result*  of  thane  works  contradict  ono  another;  la  BAny 
respect*  th*  physical  sense  of  Jho  Initial  hypotheses  la  atill  not  clear.  Thu*, 
#,gt,  Dulyko  prqoooi  that  the  atmospheric  stratification  bo  considered  within 
tho  frameworks  of  tho  simplest  system  of  th*  logarithmic  boundary  layer,  formally 
replacing  tho  Korean  "universal  constant"  by  a  variable  parameter,  a  funotion  of 
stratification.  In  Bialyko'e  system,  the  basic  characteristic  of  tho  *ubetrate, 
rougluiosa,  lo  also  a  function  of  meteorological  conditions,  Tho  purely  formal 
naturo  of  those  relations  is  ono  of  tho  shortcomings  of  Budyko's  system.  Zt 
should  alto  bo  noted  that  tho  obsorvod  profiloo  of  wind  distribution  with  height 
regularly  deviate  fro*  th*  logorltlunie  law  during  stratification  conditions 
which  dlffor  from  neutral  equilibrium, 

Lulkhtaan  proposes  s  more  complete  outhod  of  approximating  wind  and  tempore* 
turo  profllos  (an  exponential  law  with  a  variable  exponent),  which  makes  it  po*« 
sibl*  to  dlrcoru  th*  naturo  of  deviations  from  tho  logarithale  law  under  various 
condition*  of  otaospherio  stratification.  However,  Laikhtaan'*  ayotea  contain* 
too  many  fro*  paraaotoro  which  have  to  bo  determined  in  each  individual  ease. 

This  crcatos  difficulties  familiar  in  determining  these  paraaotoro  fr ca  eepirloel . 
data  and  doorcases  the  computational  accuracy. 

Those  critical  roaarka  by  no  means  are  meant  to  detract  from  tho  value  of 
tho  results  obtained  by  Budyko  and  loikhtnan  when  solving  individual  problems; 
however,  they  indicate  tho  necessity  of  developing  tho  theory  further  and  making 
tho  initial  physical  hypotheses  tore  exact, 

IThcn  analysing  tho  highly  complex  phenomena  of  ground -lay  or  turbulence, 
w.horo  tho  temperature  factors  ploy  an  cj,*...tiol  rolo,  it  is  expedient  to  uso  tho 
moth  ode  of  tho  theory  of  similltudo  which  arc  widely  used  in  applied  uurodynando* 
and  thermal  pliysics,  and  are  tho  generally  accoptod  method  of  investigation  in 
thi.i  area. 

In  1?Ij3,  A.  V.  Obukhov  attempted  to  apply  methods  of  tho  theory  of  eloill- 
tudo  to  problems  of  ground-layer  physics  £ 7 j .  Tho  results  obtained  in  this  work 
wore  subsequently  dovolopcd  by  A.  S,  lorvin  (8).  The  theory  developed  in  (?)  and 
£8]  evidently  glvo  a  satisfactory  qualitative  description  of  tho  processes. 

Furthermore,  tho  data  used  in  [7]  to  dot< rmino  the  numerical  parameters  in 
tho  propooed  systems  wore  not  sufficiently  reliablo  (tho  critical  Richardson 
number  van  ciatakcnly  assumed  to  bo  1/11,  on  tho  basis  of  Sverdrup's  data),  which 
made  It  Impossible  to  make  diroot  uso  of  tho  formulas  obtained  in  this  work  in 
actual  computations.  ’ 


Tho  proaunt  work  give#  on  analyst#  of  the  procosso#  of  turbulent  ng  l* 
tho  ground  layer  of  tho  ataosphero  on  the  basis  of  a  systematic  application  of 
tho  acthods  of  the  theory  of  similitude,  and  the  valuos  of  the  numerical  para¬ 
meters  aro  core  exactly  defined  by  using  a  sufficiently  large  amount  of  empirical 
data  on  gradient  observations,  obtalr.od  from  tho  expeditions  of  tho  U&in  Geo¬ 
physical  Observatory  and  the  Geophysical  Institute  of  the  Accdcsy  of  Sciences  of 
the  USSR.  On  this  basis,  working  formulas  were  obtained  for  computing  tho  baa la 
characteristics  of  tho  ground  layer,  vis.,  turbulent  heat  transfor,  friction, 
tho  austausch  coefficient,  and  moisture  flux,  from  gradient  measurement  data* 

The  computational  method  is  illustrated  by  specific  examples. 

1.  THE  LOG/RXTWUC  BOUKDAKS  LATER 

When  analyzing  the  processes  in  the  ground  layor  of  tho  atmosphere  on  a 
theoretical  basis,  we  will  proceod  from  the  generally  accepted  system  of  a  cur¬ 
rent  above  an  infinitely  rough  srurface  whose  horizontal  properties  are  assumed 
to  be  ruito  uniform.  The  averaged  characteristics  of  the  current  in  this  system 
are  a  function  only  of  the  vertical  coordinate  z.  The  most  important  character¬ 
istics  are  the  momentum,  heat,  and  humidity  fluxes. 

Tho  momentum  flux  can  be  treated  as  turbulent  friction  stress.  Instead  of 
turbulent  friction 

T  -  -  pu'w'"  (1) 

where  u*  and  w*  are  the  pulsations  of  the  horizontal  and  vertical  wind  velocity 
components,  p  is  air  density,  and  the  bar  indicates  averaging,  it  is  convenient 
to  examine  the  dynamic  velocity 


Within  tho  confines  of  the  ground  leyer,  t  and  the  turbulent  heat  flux  q 
con  bo  considered  to  be  practically  independent  of  height  s. 

The  condition  that  fluxos  t  and  q  are  constant  (within  the  given  toloranco) 
can  sorve  to  dotonaino  the  actual  concept  of  the  ground  layer.  Let  us  attespt 
to  give  an  approximate  estimate  of  the  height  of  the  ground  ltyer  on  the  basis 
of  possible  changes  in  T.  We  will  proceed  free  the  averaged  equations  of  hydro¬ 
mechanics  in  a  Coriolis  force  field.  The  corresponding  equation  for  ths  x-coor- 
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dinate  (Kind-volodty  direction  at  the  earth's  surfaco)  in  a  Quasi  stationary 
case  has  the  following  font 

2551  -  _  1  £2  ♦  ^7  (3) 

3s  P  a* 

whcro  3?/3x  is  tho  pressure  gradient,-  l  the  Coriolis  paramo tor  and  v  the  ooo- 
poncnt  of  averagod  wind  volocily  along  tho  y-uxis. 

Let  us  integrate  both  sides  of  tho  equation  with  respect  to  height  within 
the  licit*  of  a  layer  of  thickness  H  and  cstimato  tho  right-hand  Bidet 

d2Ual»l.  "  if .,7U<xi|iU  ft) 

P  0  p  ^  ‘  0  p  **  I 

Tho  discarding  of  tors  <v  leads  to  a  strengthening  of  the  inequality,  since 
tho  Coriolis  force  partially  compensates  tho  effect  of  tho  pressure  gradient,. 
Introducing  tho  dynamic  velocity  v#  and  tho  gcostrophic  wind  velocity  Vg  • 

-  (l/pi)  Jop/ax|,  wo  can  write  the  resultant  inequality  in  the  following  feral 

v*(0)  -  vf(H)  <Mvc  (5) 

2 

Lot  us  define  H  such  that  the  relative  change  of  v^  in  a  layer  of  thickness 
B  does  not  exceed  the  tolerance  &,  i.o., 

v*(0)  -  v*(H> 

—  v--—  £a  (6) 

v;<o> 

On  tho  strength  of  inequality  (5)  it  sufflooe  that 

cv^(O) 

w<  *  ^  •  r—\ 


in  order  that  (6)  be  fulfilled.  The  ratio  of  friction  velocity  to  geostrcphlo 
wind  velocity  can  be  estimated  to  bo  a  value  of  the  order  of  0»Q£t  , 

—  ~o.o5 
s 

freo  which  it  follows  that 


,  v 

B  C2.$  •  lO^o-f 


-  $  - 

When  Yg  10  m/see  and  4  •  Kr^socrl  wo  got 

B  ~a  •  2$0  &« 

With  •  tolorance  a  •  20%  we  get  the  estimate  of  the  height  of  the  ground  layer 
which  we  eooict 

H  •  $0  a. 

Within  the  limits  of  this  lryer,  v#  con  be  considered  practically  constant  and 
the  affect  of  the  Coriolis  force  (rotation  of  win!  with  height)  can  be  neglected. 
The  estiaato  obtained  agrees  quite  well  with  observations. 

Under  conditions  of  neutral  stratification  the  processes  of  turbulent 
nixing  in  the  ground  layer  can  be  described  by  the  system  of  the  logarithsdo 
boundary  layer.  The  corresponding  laws  have  boon  studied  in  detail  in  eaperi- 
centol .  aerodynamics,  and  are  widely  used  in  meteorology. 

Let  us  bear  in  Blind  the  derivation  of  the  logarithmic  law  of  wind  distri¬ 
bution  on  the  basis  of  the  hypotheses  of  similitude.  Let  us  assume  that  for 
values  of  t  b  t*ii  where  h^  ia  the  height  of  tho  crass  (the  characteristic  scale 
of  tho  Bicro-inhoaogeneitios  of  tho  ouhntrate),  the  statistical  characteristics 
for  relative  movements  in  a  stream  an*  *nvuriant  with  respect  to  transformations 
of  sinilitude  x'  •  kxt  y 1  ■  ky,  s’  •  kz,  t*  ■  let.  In  these  transformations  the 
half-space  z  >  0  converges,  while  the  equations  of  motion  remain  constant.  This 
factor  is  tho  theoretical  basis  for  the  accepted  hypothesis  of  siaiilitudo.  Let 
us  also  noto  that  the  natural  scale  of  velocity  v#  «  y/ x/t>  remains  invariant 
with  resect  to  the  indicated  transformations.  Let  us  examine  the  stationary 
regime  and  establish  a  ratio  of  the  difference  of  tho  averaged  velocities  at 
two  levels  zj  and  Zj  to  the  dynamic  velocity  v#.  The  corresponding  non-diaoa- 
sional  magnitude  is  a  function  of  and  Zg  and,  on  the  strength  of  the  assump¬ 
tion  of  the  self-similitude  of  the  current,  can  be  a  function  only  of  the  ratio 

▼#  ~'~T  1  • 

Lot  us  determine  the  f era  of  function  f(C).  Evidently  for  ell  three  heights 

*3  >  *2  > 


▼(*2>  -  *(Bj) 
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u(=3)  -  “(x^)  -u(Xj)  -  u(*2>  ♦  u(a2)  -  u^) 


and  along  with  this. 


*1  *2  *1 


Frca  this  it  follows  that  function  f  satisfies  the  functional  equation 

•  3^)  -  r(^)  ♦  f(^>, 

c,  -  *3A2)  (11) 

The  logarithmic  function  f(ij)  ■  C  In  £  is  the  only  solution  of  .  this  functional 
equation.  Assuming  C  ■  l/x,  wo  get 

7(t2)  -  7(0  = 

— - - —  -  -  In  (12) 

v*  •  *  *1 

whoro  x  is  the  familiar  Kansan  constant.  According  to  empirical  data,  x  asOji. 
Equation  (12)  can  bo  written  in  the  usual  differential  fora,  examining  the  in¬ 
finitely  close  values  and  ij* 

dv  Vu 


Equations  (12)  and  (13)  do  not  contain  characteristics  of  a  particular,  substrate 
but  can  pertain  to  any  substrate,  if  the  condition  Sj,  z2  t>  h1  is  fulfilled*^. 
Then,  too,  formula  (13)  dofines  only  changes  in  scan  wind  velocity  with  hoight. 
Tho  proportica  of  tho  substrate  srust  bo  considered  in  order  to  determine  the 
absolute  value  of  v(s). 

I.’ow  let  us  nosvno  that  observations  of  rind  velocity  are  conducted  at  a 
definite)  height  H  abovo  cose  definite  substrate.  Let  us  assume  that  wo  can  con¬ 
duct  independent  ocasurcaonts  of  the  turbulent  friction  and,  accordingly,  in 
ooch  individual  case  no  can  determine  vu  ■  \J -r/p  .  Tho  value  v#  can  bo  dator- 

^  Dotomination  of  the  vnluos  of  hoight  t  in  fomula  (13)  involves  a  cortain 
arbitrariness  in  tho  choico  of  the  starting  point  for  the  computation  (within 
tho  Units  of  tho  height- of  tho  gross  30*  Howover,  nhon  t  Jt  h.,  this  ini o fi¬ 
nite  ness  is  of  no  significant  value.  1  ^ 
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mined,  e.g.,  from  thcrBOancracmeter  observations  of  pulsations  u*  wd  w1,  or 
oumraril  y,  on  tho  basis  of  mersuraaent  of  tho  dr  a  2  intensity  at  the  earth's  arur- 
faco.  This  latter  method  is  used  in  practice  when  studying  turbulent  motion  la 
tubes.  Sheppard  [ 9]  attempted  to  use  the  dynamometric  method  of  measuring  X 
under  atmospheric  conditions. 

A  comparison  of  a  number  of  observations  of  v(H)  and  x  allows  us  to  doter» 
nine  the  relationship  between  these  magnitudes.  Aerodynamic  experiments  teach 
us  that  with  large  Reynolds'  masters  and  surface  "roughness"  the  dependence  of 
X  on  v  is  of  a  quad ratio  nature,  from  which  it  follows  that 

v*  -  y(H>  *(H)  (Ur) 


where  y(K)  is  a  non-dimensional  coefficient  which  is  a  function  of  the  properties 
of  the  substrate.  At  a  fixed  height  H  tho  "drag  coefficient"  y(H)  can  serve  as 
an  objective  characteristic  of  tho  properties  of  the  substrate  with  respect  to 
its  dynamic  influence  on  the  current.  However,  use  of  y(H)  has  the  disadvantage 
that  a  specific  observation  height  must  be  selected.  The  dependence  of  y(H)  on 
the  observation  height  H  can  be  easily  established  by  substituting  v(H)  •▼#/y(H) 
in  formula  (12).  For  any  two  heights  p  h^  we  will  have 


1  1 
t(h2)  ”  r^) 


1  ^2 

-  In  ~ 


h 


05) 


From  (15)  it  follows  that,  in  particular,  y(H)  decreases  with  height.  Taking 
the  or.tilogarithms  and  combining  the  magnitudes  which  contain  and  respec¬ 
tively,  ww  get 


^  (exp  -  x/y^))  •  Hg  (exp  -  VrO^))  *•  l»fl 


(16) 


i.e.,  a  magnitude  which  is  not  a  function  of  height.  Thus  the  magnitude  h^ 
which  has  length,  is  determined  only  by  tho  properties  of  the  substrate;  It  la 
called  Mynsaic  roughness."  Let  us  express  the  drag  coefficient  y(z)  by  h^* 


r(a) 


0 


(17) 


whonce  on  the  basis  of  (Ur)  we  get  the  deslrod  wind  velocity  distribution* 
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Tho  method  given  above  fcr  introducing  tho  concept  of  roughnosa  of  th®  eub- 
etrsto  has  tho  advantage  that  it  depondo  exclusively  on  tho  properties  of  th® 
current  at  rathor  great  heights,  whero  thcro  are  sufficient  grounds  for  usinj 
tho  universal  laws  of  developed  turbulcnco.  In  Dost  coses,  however,  wo  do  not 
have  at  our  disposal  tho  means  for  raking  direct  measurements  of  t  (and,  ac¬ 
cordingly,  y(H)),  ard  in  this  regard,  when  caking  practical  determinations  of 
tho  characteristics  of  dynamic  roughness,  ro  oust  use  the  properties  of  tho  wind 
profilo  which  can  bo  do  to  reined  directly  free  observations.  T7hen  dealing  with 
a  nature  vegetation  cover,  additional  difficulties  arise  in  connection  with 
choosing  tho  start  of  computation  of  z,  A  number  of  authors  (Paesehko  £10], 
Konstantinov  £11])  rccomsend  tho  uso  of  a  cortain  arbitrary  level  for  tho 
start  of  height  computations;  .this  level  lies  between  the  soil  and  tho  top  of 
tho  grass  h^.  This  level  can  bo  called  tho  hoight  of  tho  displacement  layer. 

Iho  concept  of  "displacement  height"  z^  can  bo  introduced  into  tho  general 
system  as  follows.  Equation  (13)  describes  the  asymptotic  regularities  occur¬ 
ring  when  z  ^>b^,  and  in  this  region  it  ia  insensitivo  to  slight  changes  in  th® 
starting  point  for  computation  of  z  (within  the  limits  of  the  top  of  the  grass 
hj_).  Let  us  now  oxamino  tho  region  of  values  of  z  which,  although  they  cxcood 
h^,  aro  nevertheless  comparable  with  it,  and  arc  such  that  tho  ratio  hj/x  can 
bu  treated  as  a  first-order  value.  To  bo  specific,  wo  will  compute  z  from  ground 
lovol.  In  this  case,  a  numerical  correction  factor  f(h^/x)  should  bo  introduced 
into  formula  (13);  this  describes  tho  deviation  from  the  nutcmodular  regime, 
connected  with  the  direct  effect  cf  tho  grass* 


dv 

dz 


XX 


(19)  . 


Evidently,  when  x — joo,  formula  (1?)  should  convert  into  (13),  from  which  it 
follows  that  f(0)  -  1.  Expanding  function  f  in  oorios,  w®  get 

—  y 

“  •  U  ♦  aOij/z)  -♦  ^(h^/s)2  ♦  .,,]  (20) 


Lot  us  now  introduce  a  now  starting  point  for  computations  cf  x,  assuming 
fc  -  *'  ♦  xJioro  z'  is  comparable  with  h^,  and  rewrite  tho  equation  with  re¬ 
spect  to  th®  n c.v  variable* 

^.^£i®(a-^.)  (yz’)  ♦p*(h1A«)2  ♦  ...)  (a) 


9 


Lot  ua  select  such  that  In  ajq'unrion  («'l)  tho  flrat-oNer  Un  reverts  Vo 
loro.  With  a  corrtspondii.  choico  of  with  an  accuracy  up  Vo  Vho  aeooud-onltr 
tons*,  wa  gat 


Thus,  tho  height  of  tho  displacement  lay  or  can  bo  defined  ••  tho  height  of  a  one 
arbitrary  level  of  coapuUtion,  using  v.hich  *o  4,0 1  tho  bait  approximation  of  tho 
wind  profile  by  tho  logarithaio  law  in  a  layer  situated  above  tho  grata  layer* 
Lot  ua  note  that  tho  phyaical  dotumination  given  abova  of  dynaailo  roughneta  h0 
la  insonaitiva  to  a  substitution  of  »  -  ^  for  »  (ainca  R  ^h^)|  howaver,  la  tho 
final  fcrsula  for  tho  wind  velocity  profila  no  ahould  caleulet#  tha  blight  fro* 
tho  lavol  of  tha  diaplaotaent  ltiyor,  i.e,,  replace  a  by  a  -  a^» 

v(»)  •  la  (21) 

0 


Tho  characteristics  of  tha  substrate,  a^.  can  ba  dote  rained  eapiri- 

colly  on  tho  basis  of  eoaaurcncnta  of  wind  profile  In  tho  layer  above  tha  grata 
level,  urdor  condition!  dote  to  equilibria.  To  increase  tha  computational 
accuracy  w«  should  usa  data  averaged  for  a  group  of  analogoua  oaaoa. 

Lot  ua  use,  as  an  example,  values  of  1^  and  hp  according  to  Paeschle'e  work 
[10]  (table  1)» 

Table  1 

Characteristics  of  tho  substrate 


ca 

bo,  ca 

Snow  uurfaco 

3 

o.S 

Airport 

10 

2.5 

■Sugar  boot  plantation 

6*5 

Whoat  field 

130  1 

5 

Somo  data  on  tho  question  of  choosing  the  initial  level  s^  can  ba  found  In 
on  artlclo  by  A.  R.  Konstantinov  [11],  It  is  werth  noting  that  tho  dynaalo 
roughnoas  of  a  whoat  field  Is  less  than  that  of  a  sugar  beet  plantation,  although 
the  graos  is  throo  times  highor  in  tho  first  caso.  In  the  caso  of  *  low  grass 
stand  (ctoppo)  the  value  of  does  not  play  on  essential  rule,  and  when  oco- 
puting  \  and  v#  free  observations  mode  at  heights  of  mors  than  1  meter,  ww  can 
considor  formally  that  •  0,  l.e.,  we  can  compute  tha  height  directly  fna  tha 
ground* 


-  la 


lit  Author  aeotiaiui  of  IM i  «>« W,  *>»«., u.livg  tho  offset  of  atroUft- 

e«U*n,  we  will  o<iu>lo>ur  that  height  i»  »'.0. f«v*  »i*«a  arbitrary  l«wl 
(the  Mltpleounwil  layer*),  net  I1..UW11;  11  lit  ra<y  way,  while  th§  dynaaie  WU*h* 
Data  will  be  computed  tty  mao  giv»n  uhar  ev  to  rletxO  of  the  eubstrato  wtilah  U 
independent  of  notoorelogU'al  OoiwUU.uui* 

#,  iv.sio  c(UitAOii  ni;>ria  oi‘  mi;  r,  :,w:  >nt  hwm  xn 
ak  iKiijj.oan.uius  Tcs'wuruij;  uiauu 


One  of  tho  oort  important  practice!  aU«i>adtt>i‘l»UtMt  of  l)*t  tuibuleat  ngUi 
la  tlto  grouikl  Itytr  of  iho  ebieaphera  it  llm  vertical  tuibulent  hetl  fluxi 

q  •  e^  cW?  (tli) 


whoro  «p  la  the  apwclfio  heat  of  the  u  tt  constant  pr\>»kur«,  p  la  denalty,  »• 
and  ?'  are,  respectively,  the  pulsationa  of  tho  vertical  tvuwk  velocity  eoaponaat 
*nd  of  temperature,  caused  by  U10  piui  .ago  of  turbulent  cloiaenta  through  a  cly«* 
point,  and  the  bar  indicate*  averAging,  The  nagratule  q  la  tho  average  Mount 
of  hotV  carried  by  turbulent  pulsatlona  acres*  a  unit  trot  por  unit  tint*  ITo 
have  uufflclont  ground  a  far  considering  that  for  all  intent*  and  purpoata  the 
turbulont  hoot  flux  q  in  tho  grenmi  layer  unior  atationary  condition*  is  not  a 
funotlon  of  hoight*^,  Xnatoad  of  q  wo  stay  alto  uao  tho  "tenparaturo  flux* 


sS-  -  W 


<tf> 


The  ungnitude  of  tlt>  turbulent  hoot  flux  q  «'nr.  be  determined  dlrootly  ex¬ 
perimentally,  on  tho  baoia  of  electronic  coasuruiviUa  of  tho  puloationa  of  too- 
poratur#  T1  end  of  tho  vortical  wind  velocity  component  «•»  liodom  technology 
hot  thovm  that  ouch  meaaurenento  wo  possible,  in  principle  [1?,  13] •  neverthe¬ 
less,  in  pructleo  ono  cuot  still  uao  iivlirect  method  a  to  doteraino  q,  based  on 
aix^ler  gradient  ceuaurcantta.  To  Interpret  these  tcasuroconta  correctly,  ono 
rust  tnvoatlr.ato  the  connection  between  the  characteristics  of  turbulence  q  and 

*^Rorov«o  arc  dip res sing  frox-v.  examination  of  radiational  energy  fluxoa.  Striotly 
epeoUnc,  tho  tou.1  flux  q  *  q^  is  not  a  function  of  height*  boro  q^  la  th» 

radiation  flux.  Then,  too,  in  tho  ground  layer,  ciiangcs  in' tho  radiation  flux 
q„  can  hardly  be  considered  essential,  This  quostion,  however,  should  bo  the 

oubjuct  of  apodal  invoottgatlons. 
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v#  <mi  the  distribution  of  *r,u\  wind  veK^lty  *'*1  temperature,  Whoa  solving 
this  ii'dMm  ws  «U1  follow  Uta  fe>th»t*  of  lha  tlmuiy  of  eLfttlituds  tnl 
Vo  esWltllih  «  c.-stim  oi u>  a  tilnlavn  imthor  of  pnr«;.oUiu  which  Uoiorlb*  the 
turbulent  lYgUs  in  an  irOioaogenoous  ttuperature  siudLum, 

TUs  Inhojasl^nolUes  of  Dot  U»i|>nr«Ui'«  fi  oil,  Wing  of  «  ay • toastie  nature 
(change  of  cam  temperature  with  height),  okott  *  definite  Influence  on  tho  gen¬ 
eral  turbulent  regime  (tho  offool  of  ArchlaoJoan  forces).  i’roviUud  that  tha 
twporature  pulsations  axt  alight  co&parod  with  tho  naan  t«ag>«rature  of  tho 
layer  tho  equations  for  V  dyrualoa  of  an  irhoaogenoous  temperature  medium 
oan  bo  written  la  tho  follow  „  fonsi 


Sa..l2l 

dt  P  Oat 


&  ,.i2ii 

dt  p  if 


dvr 

dt 


1 

p  do 


6u  .  dv  to  n 

da  Oy  da 


(26) 


Zn  this  system,  p^  and  indicate  deviations  from  tho  standard  values.  Tho 
etcpUficattono  mad*  when  dsrlving  tho  system  of  oquatlona  arosnoglect  of  tho 
Coriolis  for oo  aid  tho  radiation  influx  of  heat,  end  also  tho  linearization  of  tho 
standard  statistical  distribution  of  prossure  and  teeporaturo.  This  latter  in*  . 
dicatos  that  changes  in  density  duo  to  pressure  chance s  are  neglected,  end  it 
assume*  that  tho  deviations  of  density  and  teeporaturo  from  tho  standard  values 
are  proportional  (L,  0.  Landau  and  E.  Lif shits  [lkt  chapter  5])*  Thoso  sim¬ 
plifications,  used  in  tho  convection  theory,  allow  us  to  d  iscribe  tho  Archi¬ 
medean  forco  by  tho  component  (gA0)T^.  Thus,  the  equation  contains  a  diaoa- 
sional  constant  c/?0,  thich  vo  should  considor  in  the  future  when  establishing 
the  crltorla  of  similitude. 

Lot  us  no  to  that  wo  cannot  linearize  tho  equations  of  velocity  variations, 
since  in  this  case  turbulence  would  be  lost,  Zn  addition  la  tho  equations,  tho 
toraa  containing  viscosity  and  boat  cw  activity  would  bo  oalttod  %  It  It 

A  1  r  1-'~  -  -  —  -  -  —  _  -  -  - 

*  Uni  or  conditions  of  s  developed  turbulent  regime,  these  terms  must  be  eon- 


IS 


natural  to  aa;uao  tlut  churns  In  man  velocity  ant  temperature  with  holght  can 
1m  oxproaarJ  by  coonlinato  i,  panv.etu r  g/7^,  taxi  tho  "oxtonul  pa  ramotoro*  v# 
ojti  q,  while  the  corresponding  equations  cun  bo  written  In  non-dimensional  fora, 
since  they  Jo  not  oontaln  other  dimensional  constants.  This  propooition  la  the 
basio  hypothesis  of  tho  theory  or  similitude,  forsulated  in  tho  first  section  of 
the  prosent  work,  gcnoralitui  for  tho  caso  of  an  inhemogonooua  tcaporatur*  norilua. 

Tho  hypothoai*  of  tho  theory  of  similitude  which  wo  uaod  agrcee  with 
equations  (26)  and  lo  equivalent  to  tho  proposition  that  tho  system  of  equation* 
(26)  together  with  tho  condd  -Luna 

■  const. 

_  (27) 

-  puTw1"  •  x  '•  const. 


ore  an  analoguo  of  tho  boundary  conditions  and  dofino  the  statistical  character¬ 
istics  of  tho  turbulent  regime  unequivocally.  Thus,  the  throo  paramo  tors  c/T{vi 
and  q/CpP  can  bo  considered  tho  definitive  characteristics  of  tho  turbulonce 
of  tho  ground  layer  (In  tho  layer  above  the  crass).  Free  tho  so  parameters  wo 
Can  establish  unequivocally  (with  an  accuracy  of  tho  numerical  coeff icicnta)  tho 


scale  of  longth  L  and  taaperaturo  7„  which  can  be  written  in  tho  following  fonts' 

v?  ,  _ 


L  •  - 


r.  jT  * 

V  c  p 
o 


_  JL.SL 

xva  cp* 


It  ie  natural  to  uso  dynamic  velocity  v#  as  the  characteristic  velocity 
scale.  Tho  minus  sign  and  tho  Korean  constant  x  or a  introduced  for  tho  sake  of 
convcnionco.  The  signs  of  L  and  T*  arc  determined  by  the  nature  of  tho  strati¬ 
fication,  77ith  stable  stratification  tho  turbulent  heat  flux  is  directed  down¬ 


ward,  q  <  0,  and  correspondingly  L  >  0  and  >  0.  Uith  unstable  stratification, 
on  the  other  hmd,  q  >  0,  L  <  0  and  Ttt  <  0.  Thus,  wo  rust  visualize  two  quali¬ 
tatively  different  regimes,  corresponding  to  the  cases  q  <  0  and  q  >  0.  Those 
regimes  should  uni to  as  conditions  of  noutral  stratification  (q  -  0)  aro  ap¬ 


proached. 


Lot  us  examine  tho  non-dimensional  magnitudes  ~  “  and  *r-  -p  (from  now  on. 


sidcrcd  only  when  investigating  tho  vury  fine  details  of  the  microatructuro  of 
tho  wind  and  temperature  field.  Tho  vortical  transport  of  cy.mrr.tam  and  beat  is 
caused  by  tho  inhomogenoitios  of  somo  "moan  scale,"  for  which  the  direct  in¬ 
fluence  of  viscoaity  and  heat  conductivity  ore  rathor  slight. 


tho  bar  which  irvlicatos  averaging  dll  bo  omittod).  Thoso  non-dinonsionj& 
characturiotics  of  tho  ovor&cod  field  of  volocitico  and  tonporaturoe  should  bo 
definite  functions  of  tho  "external  parameters"  and  of  coordinate  t.  Tho  only 
non-dinoaaional  combination  wldch  no  can  nako  froa  q/OpP«  ▼#»  gA0  and  S  io  i/L, 
froo  which  it  follcwt  that  . 


and  • 


vr.  dv 


o^A) 


e  d? 

5"  d* 


$2(*A) 


dv 

da 


s 


c^(tA) 


(29) 

(30) 


(&•) 


(30*) 


chare  Tv  and  1  are  determined  by  formula  (28). 

Let  ua  introduce  tho  concept  of  the  austausch  coefficient.  Lot  us  aeeumo 
formally  that 


x  •  fS 


dv 

da 


<1 


V 


i  da 


(31) 


and  call  tho  dynaaic  austausch  coefficient  and  the  coefficient  of  turbulent  heat 
conductivity  K  and  Ky  respectively.  Introducing  tho  cagnitudes  v#  ■  V^x/p  and 

T*  ■  -  — —  -2—  in  place  of  x  and  q,  and  using  equations  (29)  and  (30),  we  got 

*v#  °pp 


E"  TOTT*  *Tm5p&Z7 


02) 


J.’o.t  lot  us  exaaino  tho  hypothosts,  aharod  by  a  najoxlty  of  meteorologists,  that 
within  tho  limits  of  Doteorological  observations  wo  can  consider  that  K  •  , 

**  Conorally  speaking,  K  >  since  tho  offoct  of  proscuro  pulsations,  as  well 
as  BlxinCf  can  bo  expressed  in  a  acao&tua  exchange.  IJowovor,  as  of  new  wo  have 


fra i  which  it  follow*  that 


<^(zA)  “  <g(=A)  *  «\ *A) 


(33) 


Tho  similitude  of  tho  tcnporaturo  and  wind  profiles  follow*  dirootly  froa  th* 
accepted  hypothec!*  that  X  -  Kj.  Dividing  (30)  by  (2?)  w*  get 


dT 

dv 


xT_ 


V 


<310 


and,  accordingly,  for  any  heights  and  Hg 


T(H2)  -  TfH^)  -  ~  Rf^)  -  7(5^)) 


(35) 


Thus,  tho  ratio  of  the  difference  of  mean  temperatures  at  two  levels  and  Eg 
to  the  difference  in  velocities  at  the  samo  heights  is  not  a  function  of  choice 
of  heights  and  Hg,  but  is  determined  entirely  by  external  conditions  -  th* 
ratio  of  tho  turbulent  heat  flux  q  to  turbulent  drag  resistance  t. 

let  us  nor  shenv  that  tho  non-dimensional  factor  c(zA),  where  L  • 


v? 


ft  O 

X  «*■ 

*0  V 


,  is  directly  connected  with  the  Richardson  number  at  a  given  level. 


Substituting  tho  values  dv/dz  and  dT/dz,  determined  from  formulas  (29)  and  (30), 
in  the  expression  for  the  Richardson  master 


Ri 


d? 

p  c: 
(dv/dz)2 


(36) 


VO  got 


no  conv.e'ing  ovidenco  that  thio  difference  is  essential.  Tho  theory  developed 
in  tho  present  work  can  bo  generalised  fer  tno  case  iw’a  *  a  f  \  if  w*  jnoplsc* 
T  by  T/u  in  all  instances. 

It  follows  that  T  should  indicate  potential  te.tpu  nature,  sine*  T  does  not 
chango  with  vertical  ohifts  of  the  turbulent  olcwtts  (the  stuto  of  tho  latter 
can  bo  considered  adiabatic).  In  tho  ground  layer  tho  mr.cncal  values  of  po¬ 
tential  and  molecular  temperature  are  very  close.  ”.’ith  tho  largo  tenperatur* 
gradients  usually  observed  in  tho  ground  layer,  tho  difference  between  tho  gra¬ 
dients  of  potential  and  mol ocular  temperature  are  inconsequential;  however,  la 
•tates  cloo*  to  ioothercy,  this  difference  is  significant. 


-  is  - 


Ri 


**»8 


0  Vw9(sA) 


or,  using  the  determination  of  the  scale  of  L  (28), 


Ri 


a  1 

L  *  vCi/L y 


(37) 


(30) 


fr oa  which  it  f  oUero  that  the  dcpcndonco  of  the  Richardson  nuaber  on  height  is 
defined  by  &  single  parameter  -  tho  scale  of  L. 

Ccapering  formula  (32)  for  tho  austausch  coefficient  with  the  expression 
for  tho  Rich  aid  son  nuaher,  wo  got  on  important  relationship  bou.’een  the  auetausdh 
coefficient,  the  scale  of  L  and  tho  Richardson  nuahert 


K  -  xv#L  .  Ri 


09) 


Lot  us  explain  the  physical  sense  of  "the  scale  of  L."  Under  any  con¬ 
ditions  of  stratification  we  have 

—  <?(*/l)  (i*o) 

da  xs 

Let  us  fix  value  z  and  decrease  aagnitudo  q  infinitely,  approaching  tho 
conditions  of  neutral  stratification,  vrhich  corresponds  to  infinite  growth  of' 
the  scslo  of  L  (with  respect  to  absolute  nagnitude).  Obviously,  within  tnle 
range,  we  should  obtain  formula  (22),  froa  which  it  follows  that 

9(0)  -  1 

Under  given  external  conditions  characterized  by  cagnitudos  v,  and  q  an d 
tho  corresponding  magnitude  of  L,  in  tho  region  of  values  of  z  which  are  quite 
email  co rp axed  to  L,  o(zA)  will  be  quite  closo  to  unity.  This  indicates  that 
austousch  conditions  with  i  diffor  little  froa  austausch  conditions  in  « 
noutrally  stratified  ataosphore  and,  accordingly,  turbulence  is  caused  mainly 
by  puroly  dynamic  factors.  Thus,  the  scale  of  L,  first  introduced  by  Obukhov 
[7],  is  oa  important  phyoictl  charactcrioUo  of  tho  state  of  the  ground  ley  or 
and  can  ho  called  the  bedsit  of  the  substrate  of  dynamic  turbulence.  On  the 
strength  of  the  fact  that  9(0)  •  1  and  formula  (38),  when  1— *0,  we  get 


16  - 


i  - 
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This  formula  can  servo  as  the  basis  Tor  determining  tho  scale  of  L  from  empiri¬ 
cal  data  (free  the  wind  and  tesperaturo  profiles). 

Tho  function  9(2 A)  should,  in  the  general  case,  bo  determined  freo  the 
aggregate  of  empirical  data.  It  should  be  noted  that  tho  data  available  at 
present  aro  insufficient  to  determine  function  9  roliably  in  a  sufficiently  wide 
range  of  changes  of  the  argument  z/L.  However,  a  number  of  important  problems 
can  be  solved  for  the  case  z/L  <  1,  whore  wo  car.  limit  oursolvea  to  the  first 
tens  of  function  9  expanded  in  sorios.  This  case  requires  special  examination. 


3.  DET3U.2HATT  ON'  OF  TICS  TURBULENCE  CHARACTERISTICS  FRCU 
DATA  ON  GRADIENT  hEASURELENTS 

In  the  case  J  z/L  ]  <  1  we  can  limit  oursolvcs  to  the  first  terms  of  the 
function  9(2/1)  expanded  in  *  power  series  and  assume 

9(2/0- !♦? J  ,  02) 


where  ?  is  some  universal  constant  which  can  be  determined  on  tho  basis  of  em¬ 
pirical  data.  From  formulas  (29),  (30)  and  (li2),  integrating  with  respect  to  *, 
we  got 


v( 2)  -  —  [In  (z/h  )  *  p(sA)J 
x  0 

T(2)  -  T(ho)  -  [In  (zAo)  ♦  £(iA)) 


(1*3) 


Kero  wo  replaced  tho  components  p[(s-h0)AJ  by  ?(z/L),  with  the  intention  of 
using  formula  (li3)  only  whoa  z  £>  he. 

Let  us  note  that  analogous  formulas  can  be  used  to  describe  ihc  profiles  of 
tho  concentration  of  any  passive  substance  in  the  ground  lavor  of  the  atmosphere. 
For  example,  with  a  stationary  turbulent  reg ico  and  tho  absence  of  phase  con¬ 
versions  of  humidity  in  the  atmosphero,  tho  vertical  moisture  flux  ("rate  of 
evaporation")  £  -  pw'Q*‘  (Q  is  specific  humidity)  can  bo  considered  independent 
of  height  and,  analogously  to  (30),  wo  can  sat 


dQ 

da 


Q# 

»  9(»A)« 


Q* 


1  c 

*v^  P 


(U*> 
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<5(0  -  Q(h0)  -  Q(la  (z/ho)  ♦  p 


Finally,  the  expression  Tor  the  . austausch  coefficient  X  -  kv«L  •  Ri  in  approxi¬ 
mation  (42),  according  to  formula  (38),  will  have  the  fora 

SU>  -  — r  '  (W) 

1^£- 

tfith  equilibria  stratification  (|Lj»  co)  wo  got,  frco  (43),  the  usual  log¬ 
arithmic  foraulos  for  wind  and  temperature  distribution  with  height.  Non-equi¬ 
librium  of  stratification  is  described  in  (43)  by  the  component  p(z/[.)  and  leads 
to  regular  distortions  of  the  logarithmic  law.  YJith  unstable  stratification 
(L  <  0),  intense  turbulent  nixing  leads  to  equalization  of  wind  velocity  In 
different  layers  of  the  atnasphcro,  so  that  the  wind  velocity  should  Increase 
with  height  core  slowly  than  in  the  case  of  neutral  stratification,  i.e., 
should  bo  loss  than  soro.  •  Accordingly,  0  >  0.  ^ 

Formulas  (43)  for  v(z)  and  T(z)  aro  in  good  quantita-i* 

Stability 

tlvo  (and,  with  correct  selection  of  the  parameters,  also  <  / 

qualitative)  agreement  with  the  observed  profiles  of  wind  J  yf 

velocity  and  temperature  in  the  ground  layer.  Actual  nea- 
curemonts  confirm  the  presence  and  nature  of  regular  do-  ^^Inataolllty 
viations  of  the  logarithmic  law  in  the  wind  and  tempera-  *1*4’ 
ture  distribution  with  height,  indicated  by  formulas  (43).%*  Stability 
This  can  bo  seen,  o.g.,  from  the  data  of  table  2,  which  I.  /  S 

shews  wind  profiles  averagod  by  groups  with  an  identical  I  /  / 


'X 


Instaolllty 


stability  parameter  S  •  -y  (uaken  from  data  of  tho  Ualn  If 

V  A.  W 

Geophysical  Observatory  expeditions  of  1945  [15],  1947  ’S.InatabiUty 

[16]  and  1950  (17]  and  tho  expedition  of  the  Geophysical 
Institute  cf  tho  Academy  of  Sciences  of  tbo  USSR  in  19>1  J. 

(16]).  Tho  form  of  profiles  v(s)  and  T(t),  in  agreement 

with  formulas  (43),  is  d-von  in  figoro  1.  Figures  2  and  2?Sri*airf*tc5«£ 
3  Civ®  tho  avoraged  profilos  of  wind  velocity  and  tea-  ture  profiles, 

poraturo  obtalnod  by  tho  1951  expedition  of  tho  Geophysical  Institute  of  the 
Academy  of  Sciences  of  tho  USSR*^. 

The  straight  dashod  linos  in  figures  1,  2  and  3  correspond  to  the  logarithmic 
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Table  2  (conclusion) 
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Halstead  [I?]  proposed  that  tho  influonco  o f  stratification  bo  cccputod  by 
introducing  correction  factors  into  tho  logarithmic  fonauloo,  analogous  to  (1*3)* 

but  without  analyzing  tho  coofficiento  frea 
tho  point  of  view  of  tho  theory  of  eiaili- 
_ _ 'i— I,  tuda.  £ 


- -  £J«  t»  I*  ~  X.UJJ 

» 

Figure  2.  Averaged  wind  profiles  ac¬ 
cording  to  capirical  data.  • 


Figure  3.  Averaged  temperature 
profiles  according  to  empirical 
data. 


Approximating  tho  measured  wind  and  temperature  profiles  by  formulas  (ia3 ), 
wo  can  doterair.o  tho  t'urbulcnco  characteristics  froa  gradient  measurement  data. 
In  practice,  during  such  an  approximation  we  oust  first  determine  tho  starting 
point  for  coeluting  height  -  tho  thickness  of  the  displacement  layer.  Tho 
magnitude  can  bo  dotcmlnod  experimentally,  eo  that  on  tho  groph  with  tho 
loguTithsic  ocolo  tho  wind  profiles,  corresponding  to  eases  of  equilibrium  at rat¬ 
ification  (i.o.,  actually,  to  cases  of  isothersy)  would  be  depicted  by  straight 
linoo  with  respect  to  height.  Extrapolating  the  obtainod  roctilincar  wind  pro- 
Xilo  craph  to  zero  velocity*  *°  bho  value  of  the  roughness  height  hg. 

Tho  cagnitudo  tr<4  tho  parameters  v^/x  and  0/L  which  ontor  into  fomulas 
(1*3)  can  bo  cost  Accurately  determined  by  using  tho  loast-cquaros  method  to  pro- 
coss  tho  wind  profiles  measured  at  tho  saso  station,  gonerally  sparking,  under 


various  cor*ii  lions  of  stratification, 
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ThuA.  eiiUU.ing 


V1(#)  -  Al(y  *  log  *)  ♦ 


for  oach  profile,  where  1  la  tho  nuibor  of  the  profile,  ar\U  uelooting  A^,  y# 
aid  C<  because  of  tho  requirement  that  the  ouu  of  tho  squares  of  tho  deviations 

a* 

bo  cinioua. 


Z  I*i(T  *  1«S  *t>  •  C1>k  -  WJ* 


0*7) 


wo  got  for  each  profile 


K 


In  10 


£  -  -1  la  10 

U  *4 


and  wo  got  a  cocaon  roughness  height  b0  *  10**1*  for  all  profiled* 

Having  dc  to  reined  p/L  for  each  profile  by  tho  indicated  nwthod,  knowing  h^, 

.  «5>  •  «3> 

£  tho 


and  ccoputlng  tho  value  of  tho  stability  parameter  S 


w» 
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which  follows  froa  (ii3).  Tho  nuabor  p  con  bo  detcinincd  as  tho  r ogrosslon 
efficient  of  values  of  F(p/L)#  coraputod  froa  the  pi  ovioualy  oalculatod  fi/lt  for 
the  computed  valuoa  of  S,  Tho  regression  coefficient  P,  coaputed  froo  tho  data 
of  the  four  expeditions  listed  in  tablo  2,  is  0,62;  the  accurocy  in  dotoraining 
P  in  this  coso  is  probably  not  noro  than  102.  A  determination  of  fi  from  tho 
data  of  Just  ono  L'ain  Geophysical  Observatory  expedition  [16J  yielded  a  value 
of  0.57. 

Using  formulas  (li3)  wo  can  eoeputo  tho  drag  velocity  v#J  as  well  as  tho 
turbulonco  characteristic  which  has  tho  cost  practical  interest,  i.o.,  tho  heat 
flux  q,  using  tho  results  of  wind  velocity  and  tenperaturo  coasuraaento  at  only 
two  heights.  For  exncpla,  lot  z^  -  K/2,  12  *  H(  and  *jj  •  2H,  and  lot  us  as ruas 
that  tho  valuoa  Tj  »  1(6^),  7j  -  T(sj),  and  Vj  •  v(*2)  »/«*o  have  boon  measured. 
Tbon  from  (1x3 )  wo  got 


-  r?  - 


y  «•  «—«* 
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V  ^ 
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.  ,«v  (T,  --,)  ,  ...  Yf<T>  -  V  ~i_ 
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1  «  O.fiVji  oat). 


litre  wo  uaed  tha  value  »  ■  Q.VJ  for  the  K.>n.*n  constant,  It*  ekcnituda  H/l  la 
deten&lned  f rua  relationship  (U),  which  eatucua  the  foru 

M  A  **•»/<>.«♦  «« 


lo;  ** 


T.  -  V, 


9  *1 

whew  U  •  0,10?  II  (log  a^/ll)  •  >*-;»•  -  ,  whila  It  Ukl  V|  »ra  expreaaod  1a  Mtin 


orU  a/uco,  ivhpoctivuly. 

The  influence  of  atreti fitaUon  on  the  laajnitude*  b#  artf  q  la  exproaaod  by 
the  eppcartnco  of  tha  coeponcnta  with  It/I.  in  t ha  dcnoainaiora  of  fonaulaa  (U$). 
Aa  k  rule,  tho  corroeUon  for  ^ratification  api-twa  to  ba  alicht  (HA^IO”*), 
which  la  natural,  air.co  turUilcnce  in  the  lu'ier  part  of  tha  tround  layer  la  da* 
tenained  udnly  by  dyiuuaio  f-vUre, 

Feruulaa  (i,?)  u«i  (liO)  c-u  bo  umful  when  t~»#a-procoaain£  tbo  eradieot 
uoaour  Jtont  data.  In  bjooHic  caaca,  v*th  ha  ani  H,  theao  foraila*  ac¬ 

quire  a  relatively  aiaplu  form.  For  vxauplo,  nh>-n  ^  «  1  ca  and  U  -  1  a,  ere  have 

m  m 

B  .  q.u  <  i  .  -o.X5  «  ~  (I  ♦  /TTVX5)  •* 

v;  & 


°.09S  vn  v  (?.  -  T.)  . 
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Exanploa  of  cooputnUorj  of  tho  turbulent  hoot  flux  q  (froa  data  of  tha 
Ooophyuicftl  Institute  of  tha  ;cadu=y  of  Scioncoa  of  the  US-SJt  expedition  of  19£l} 
are  civon  in  [IB).  CoapuUtiona  using  apaoific  data  a  ho*  that  tha  aoala  of  l  la 
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u*v,\l)y  of  tho  oi-dor  of  10  o,  and  It  approaches  3— 1»  a  only  in  opeoiflo  cotta 
tilth  gmat  Instability  or  abrupt  invoraiona.  In  caaea  oloao  to  laothoray,  L 
rcachca  value*  of  several  ten*  of  actors.  Tho  drug  velocity  v#  la  About  8£  of 
tho  uiiul  velocity  at  8  n  with  uni  table  stratification,  and  about  $%  with  at  able 
atratif leution.  In  Keiakhatan,  In  eusaor,  the  turbulent  heat  flux  q  reaohao 
0,;'$  -  0,35  c*l/coMaia  on  hot  sunny  days,  while  it  la  of  tho  order  of  0.06  oal/ 
ca?tin  at  right. 

Considering  that  aoao  researchers  uto  tho  formulae  propoaod  by  Bv  ,-ko  [3) 
and  Lalkhu&n  (5)  when  determining  the  turbulence  characteristics  from  gradient 
measurement  data,  let  ua  dorivo  the  relationship  between  the  acale  of  L  and  the 
bade  par  abater  a  of  the  Budyko  end  Laikhtaaa  formulas.  Budyko  approximates  the 
wind  profile  a  by  the  logarlthado  laws 


v(i)  - -lnra 


(52) 


where  a  la  a  paraaetor  which  la  u  function  of  atmospheric  stratification  (with 
neutral  stratification,  a  rave  .a  to  unity).  Equating  tho  expreaeleaa  for  vCiji/ 
v(a^),  eoaputed  frea  fomulaa  (13)  and  (52),  wo  get  the  ratio 
£  In  u  •  la  (s^/t^) 


L  \  u  r  -  *8  ^  "  «2>  11 
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(S3) 


failing  to  the  Halt  whon  *j — >  *i  *  A  we  get 

In  b 


>1 


1  -  In  o  ♦  In 
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Lai kh tain  approximates  the  wind  profiles  by  tho  exponential  lain 

a*  -  h? 


v(t)  •  vft^)  ~ 


-  b' 


($1*) 


(SS) 


*1  -o 

whore  6  is  a  paraaetor  which  la  a  function  of  atmospheric  stratification  (with 
neutral  stra  till  cation,  b  reverts  to  tero).  Equating  tho  expressions  for  rftj)/ 
v(a^),  computed  from  fomulaa  (L2)  and  (55),  wa  get  the  relationship 
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(aj.bj)  la^-  (nj*bj)  la^. 
*2  <*l-bo°>-*l  <“£-*# 


(56) 


Passing  to  the  Holt  when  «  K,  »e  got 


6  In  rf  -  [1  -  (h  A)6)  ' 

o _ _ 

1  -  6  -  (hg/U)6  . 


(57) 


Taking  advantage  of  the  fact  that  tho  value  of  6  is  insignificantly  snail,  and 
expanding  tho  rlghthand  side  of  (57)  in  ocrics  according  to  tho  6-exp  on  on  t», 
w a  gat  the  approximation 


5 _ * 


2*0 
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(58) 
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Pros  formulas  (30)  and  (1*0)  it  follows  that  in  a  stationary  turbulent 
ground  lay  or,  the  wind  and  temperature  profiles  can  be  described  using  one  uni¬ 
versal  function  of  z/L.  Thus,  integrating  (30)  with  respect  to  z  and  sotting 


*(0  -  ; 


ste?j££ 


,  wo  get 


v(z)  -  (f(zA)  -  f(  yDj 


(5?) 


T(s)  -  T(ho)  ♦  7#  Cf(zA)  -  f(h0A)) 


In  the  preserit  section  wo  will  investigate  the  form  of  the  universal  fac¬ 
tion  f(zA)  taken  as  a  whole. 

Sinco  oCt) — yl  when  K—*0,  with  small  zA  the  function  f(>A)  io  of  on 
asymptotically  logarithmic  natuR.-t 


f(zA)  ~ln 


♦  const. 


when  ~  1 

Is 


(80) 


Hith  largo  z/l>  tho  asymptotic  behavior  of  function  f(zA)  will  differ  in  coses 
of  unstable  (L  <  0)  or  stable  (L  >  0)  stratification,  since  is  these  cases  tbert 
are  actually  t.*o  qualitatively  different  regimes  of  turbulent  motions. 

To  ealyzo  the  case  of  unstable  stratification,  first  lot  us  examine  the 
limiting  case  of  purely  thermal  turbulence  (with  no  wind).  In  this  case,  due  to 
the  lack  of  on  averaged  wind,  the  friction  stress,  on  an  average,  will  be  zero 


(v#  •  0),  while  the  turbulence  regimo  is  characterised  by  only  the  parameters 
q  and  g/T0  (the  turtulonco  receives  its  enorgy  exclusively  from  the  instability 
energy,  and  therefore  is  a  function  only  of  the  degree  of  instability,  char* 
octerized  by  tho  heat  flux  q  >  0  and  of  the  magnitude  of  the  Archimedean  force*, 
characterized  by  the  parameter  c/*0)« 

Tie  cannot  fora  the  scale  of  lencth  fron  the  parameters  q  and  g/Tej  there¬ 
fore,  tho  regime  of  purely  thermal  turbulence  is  autonodular  ,  i.*., 
its  characteristics  are  combinations  of  q,  g/T<>j  wi  t.  Proa  the  concept*  of 
■  dimensions  we  got 

T(a)  -  ♦  -far  (q/cpp)2/3  (gx/l^)"1^3  (fil) 

_t/3 

where  C  is  the  non-din ensional  (universal)  constant,  the  factor  x^  ie  intro¬ 
duced  for  convenience,  and  T^  is  a  constant  which  has  a  temperature  dimension. 

Fras  (61)  it  is  evident  that  v/ith  an  increase  in  height  the  distribution  of 
temperature  approaches  isotheTy*^.  This  is  natural,  since  in  the  case  of  un¬ 
stable  stratification  at  grea'.  heights,  large  turbulent  elements  develop  (whoso 
dicensions  are  limited  only  by  the  distance  to  the  earth's  surface),  bringing 
about  voxy  intense  mixing  of  the  air,  which  leads  to  an  equalization  of  the  tem¬ 
perature  profile. 

From  (61)  it  follows  that  the  austausch  coefficient 


*-  -=%?•  \  <<i/V>V3  <b/vV3  <x*)U/3  (“) 

V5I 

rapidly  increases  with  height,  which  is  explained  by  tho  augmentation  of  the 
turbulent  elements  with  an  increase  in  height  and  the  simultaneous  increase  in 

up  \ 

the  intensity  of  the  pulsations  • 

Formally,  formula  (61)  con  be  written 


T(s)  -  7(ho) 


C  (z tvr™  -  c  (b t/L)’ 


'  In  (61)  wo  aro  speaking  of  tho  approach  to  "potential  isothexmy"  with  an  In¬ 
crease  in  height  (see  footnoto  on  page  lU). 

y7ho  concepts  here  presented  on  the  regime  of  purely  thermal  turbulence  agree 
with  the  cyslcs  proposed  by  A.  A.  Skvortsov  [JO],  with  tho  sole  difference  that 
Skvortsov  introduces  a  concept  of  tho  discrete  cpoctrum-of  the  scales  of  tur¬ 
bulent  formations,  while  in  the  system  presented  here,  the  spectrum  of  the 
scalos  le  assumed  to  be  continuous. 
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oo  that  in  tho  caso  of  purely  thermal  turbulence,  the  universal  function 
(determined  with  as  accuracy  of  the  constant  component)  has  tho  fora  f(sA)  m 
-  C(s ♦  const. 

The  caso  cf  purely  thcraal  turbulcnco  can  be  derived  from  the  general  case 
of  unstable-  stratification  by  passage  to  tho  Unit  with  v#  — >0.  Here  L  — > 0  and. 
i/b— a- oo.  Therefcro  tho  asymptotic  behavior  of  the  universal  function  f(cA) 
is  determined  by  tho  relationship 

f(sA)~  C(zA)”1^3  ♦  const,  when  7  $  -1.  (6i*) 

I# 

This  result  indicato9  that  at  great  heights  2  £ jL]  (in  tho  ground  layer)  the 
turbulent  regimo,  in  the  case  of  unstable  stratification,  i3  dete rained  nainly 
by  thcraal  factors  (the  wind  profile  is  smoothed,  and  turbulence  receives  its 
energy  nainly  froa  tho  energy  of  turbulent  instability,  not  fron  tho  energy  of 
average  notion). 

An  explanation  of  tho  asymptotic  behavior  of  the  function  f(z A)  when  t 
in  tho  caso  of  stable  stratif  ication,  requires  th.tt  wo  ir.troduco  additional  con¬ 
cepts.  Turbulence  degenerates  in  the  limiting  case  A  abrupt  inversion  with  a 
vanishingly  vcok  wind.  The  existence  of  large  turbulent  elements  bocomce  ia- 
possible  in  tho  case  of  stable  stratification  (since  they  must  expend  too  much 
energy  on  opposing  the  Archimedean  forces),  and  turbulence  can  exist  only  in  the 
fora  of  saall  eddios.  Large  waves  cannot  lose  stability,  which  is  natural  frea 
the  point  of  view  of  the  theory  of  stability.  In  this  caso  turbulent  exchange 
between  different  atmospheric  layers  as  hampered  and  turbulence  takes  on  a  local 
character;  at  rather  high  altitudes  z  L  (or,  to  put  it  cnother  way,  with 
strong  stability,  i.e.,  at  low  heights  L  >  0)  tho  turbulcnco  characteristics 
evidently  cannot  be  functions  of  the  distance  1  to  the  substrate.  This  pertains, 
in  particular,  to  tho  nixing  coefficient  K  and,  accordingly,  also  to  the 
Richardson  number  Hi. 

Thus,  wo  cay  considor  that  in  tho  case  of  stable  stratification  with  an  in¬ 
crease  in  height  z  (or,  with  an  increase  in  stability,  i.e.,  a  doorcase  in  L), 
tho  coefficient  of  nixing  K  and  the  Richardson  number  Ri  tend  toward  certain 
constant  values.  This  is  natural,  since  with  an  increase  in  stability,  K  evi¬ 
dently  cannot  increase,  whilo  Ri  cannot  docrcaso.  Accordingly,  there  is  a  (uni¬ 
versal)  value  R  of  the  Richardson  number,  which  is  such  that  when  %/L  £  1, 

•  Ri  ‘‘-R  ■  conat,  X  -*xv#L  •  R  (6£) 
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Tho  limiting  value  or  It  evidently  cannot  bo  creator  than  tho  critical  value 
JUcr,  but  since,  asymptotically,  X  ^  0,  i.e.,  turbulcnco  does  not  completely  de- 
genorato,  R  should  bo  loss  than  Ricr.  The  limiting  value  obtained  will  be 
called  tho  stationary  Richardson  cumber. 

From  (£5)  it  follows  that  when  s/L  ^  1,  f'(0  should  approximate  l/ft,  or 

•  f(sA)  7  ♦  const.  (66) 

A  X# 

Haro  we  have  ' 

v(*)  *—  -  j  =?*  “^r  ♦  const.  (67) 

KroVv‘ 

T(s)  ~  ~  (<l/cp)J  -T-  ♦  const,  (68) 

R  o  p  v£ 

Our  formulas  (60),  (fli)  and  (66)  show  the  behavior  of  function  f(C)  when 

1,  6^-1  and  K  ^1,  respectively. 

For  an  empirical  determination  of  the  universal  function  f({)  in  a  suf¬ 
ficiently  broad  .range  of  changes  in  the  parameter  6,  us  in;  the  data  of  the  four 
expeditions,  given  in  table  2,  -and  determining  v#  and  L  (when  p  ■  0,6)  for  each 
wind  profile,  we  construct  the  empirical  universal  function 

—  !▼(*)  -  v(|ll/2)J  -  f(*A)  -  f(ll/2) 

where  the  plus  sign  corresponds  to  stable  strati fioation,  and  the  minus  alga  to 
unstable  stratification. 

Tho  empirical  points  obtained  are  plotted  on  the  graph  in  figure  h.  The 
graph  gives  convincing  evidence  of  the  suitability  of  the  hypotheses  of  simili¬ 
tude  used  in  the  present  work;  these  hypotheses  reduco  to  the  existence  of  a 
single  universal  function  f(z/L)>  Tho  empirical  points  lie  along  saooth  curves 
with  a  very  small  scatter,  despito  the  inaccuracies  of  the  wind  measurements  sad 
tho  corputation  of  L  and  v#  by  the  approximation  methods  shown  above.  Some 
scattor  of  the  points  is  noted  only  in  highly  stable  cases.  Tho  drawing  shows 
the  limiting  behavior  of  the  curve  quite  well  for  the  case  cf  high  stability 
(approaching  a  linear  profile)  and  high  instability  (aporoaching  a  constant).  . 
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